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ABSTRACT 

Stellar population and the interstellar gas-dust medium in the vicinity of the 
open star cluster vdB 130 are analysed using optical observations taken with the 6-m 
telescope of the SAO RAS and the 125-cm telescope of the SAI MSU along with the 
data of Spitzer and Herschel. Based on proper motions and BV and JHKs 2MASS 
photometric data, we select additional 36 stars as probable cluster members. Some 
stars in vdB 130 are classified as B stars. Our estimates of minimum colour excess, 
apparent distance modulus and the distance are consistent with young age (from 5 to 
10 Myrs) of the cluster vdB 130. We suppose the large deviations from the conventional 
extinction law in the cluster direction, with Ry ^4 — 5. The cluster vdB 130 appears 
to be physically related to the supershell around Cyg OBI, a cometary CO cloud, 
ionized gas, and regions of infrared emission. There are a few regions of bright mid- 
infrared emission in the vicinity of vdB 130. The largest of them is also visible on Ha 
and [Sll] emission maps. We suggest that the infrared blobs that coincide in projection 
with the head of the molecular cloud are Hll regions, excited by the cluster B-stars. 
Some signatures of a shock front are identified between these IR-bright regions. 

Key words: ISM: kinematics and dynamics - ISM: clouds - ISM: lines and bands - 
ISM: Hll regions - infrared: ISM - open clusters and associations: individual: vdB 130 


1 INTRODUCTION 

Star formation picture drawn by modern observations is 
much more intricate than it used to be a few decades ago. 
Our current view of the early phases of pre-stellar and stel¬ 
lar evolution has been expanded significantly by Spitzer, 
Herschel, HST, Chandra, XMM and other space-based and 
ground-based facilities. 

Now we have a much better understanding of sequen¬ 
tial star formation triggered by turbulence, feedback from 
outflows, supernova explosions and expanding Hil regions. 
These and other factors work together sculpting a fine struc¬ 
ture of molecular clouds with filaments, pillars, blobs, pep¬ 
pered with young stellar objects (YSOs). Given this com¬ 
plexity, along with large-scale studies, a detailed look at 
some specihc objects and regions can be useful for distin- 
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guishing between general and particular features of the star 
formation process. 

In this paper we present an investigation of the stellar 
population and gas-dust medium in the region located in the 
north-west border of a supershell surrounding the Cyg OBI 
association (Fig. [1^). 

The Cyg OBI supershell with an angular size of 3° x 4° 
is a part of Cygnus X region, the richest a nd most massive 
complex of star f ormation in the Galaxy. I Schneider et al.l 
(| 20061 . 12 OO 7 I . l201lh showed that the Cygnus X South region 
contains a large scale network of molecular clouds with ubiq¬ 
uitous recent and ongoing star formation regions. The area 
includes several OB associations and stellar clusters. Their 
energy injection from UV radiation, stellar winds and SN 
explosions could be a driving source for the supershell for¬ 
mation and expansion. The cloud complex at the south-west 
of Cygnus X has an extended shell-like structure, which co¬ 
incides well with the infrared, optical, and radio supershell 
encircling Cyg OBI and Cyg OB9 associations. Study of the 
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Figure 1. (a) A north-western part of the supershell surrounding the Cyg OBI association as it appears in t he op tical on the E-DSS2 
map. Black rectangles show the cluster region (top) and the region studied in the work of I Arkhipova et ahl (|2013l ) (bottom). A white 
rectangle indicates the region for which the PV diagram is constructed (see Section 4.1). Numbers mark Cyg OBI stars with strong 
stellar winds: 07e HIP100173 (1), 08 HD193595 (2), and Of HD 228841 (3). (b) The same region as seen in the infrared on the 3.6 Atm 
Spitzer image. Overlaid are contours of ^^CO J=l-0 emission. The box limits the area with CO-isocontours from I Schneider et al.l (l2007lV 


triggered star formation within walls of this shell can help to 
reveal the spatially resolved star formation history in the en¬ 
tire region. In this series of papers we plan to present results 
on the stellar content, young stellar objects, and the diffuse 
matter in the area, putting them into a general context of 
the star formation in giant molecular cloud complexes. 

In this work we analyze a small region in the north¬ 
western wall of Cyg OBI supershell. It contains the 
young open cluster vdB 130 characterized by an anoma¬ 
lously strong and non- u niform reddening Jr acin 3 1 19741: 
Moffat Schmidt-Kalerl Il976l : iLe Duigou Knodlsedeil 
2 OO 2 I ). the reflection nebu la vdB 130 of the same 
name (Ivan den BerghI 19661 ). also referred to as GN 
20.16.0 (jMagakianl 12003 ). multiple compact and extended 
infrared (IR) sources, and a molecular clou d of a comet ary 
shap e seen in CO emission (cloud A in I Schneider et al.l 

12007IV 

The goal of this study is to revise the stellar content 
of the cluster vdB 130 (in a view of new optical and IR 
data) and to analyse the physical relationship between vari¬ 
ous components of stellar, gaseous, and dust populations in 
the region of sequential star formation. 

This paper represents a follow up of our previous studies 
of triggered star formation regions in the n orth-west wall of 
the C yg OBI supershell. In the work of lArkhipova et al.l 
1 I 2013 I) a family of cometary globules with a few patches of 
star formation has been investigated in a region, shown with 
a lower black rectangle in Fig. [1^. 

Section 2 describes the spectroscopic, photometric, and 
interferometric observations. Section 3 presents the results 


of our revision of the stellar composition of vdB 130 based 
on the data of the 2MASS survey, the UCAC4 catalogue of 
proper motions, and the spectral types that we determined 
for a number of stars. In Sections 4.1 and 4.2 we analyse 
the structure, kinematics, and the emission spectrum of the 
interstellar medium (ISM) in the vicinity of vdB 130, and 
possible links between various components of the stellar pop¬ 
ulation and the ISM in the region. In Section 4.3 we discuss 
the nature of the IR radiation of the interstellar medium in 
the neighbourhood of vdB 130 based on Spitzer and Her- 
schel archival data. The concluding section summarizes the 
main results of the study. 


2 OBSERVATIONS AND DATA REDUCTION 

This study is based on optical observations performed with 
the 6-m telescope of Special Astrophysical Observatory of 
Russian Academy of Sciences (SAO RAS) and the 125- 
cm telescope of the Crimean Laboratory of Sternberg As¬ 
tronomical Institute, Lomonosov Moscow State University 
(SAI MSU) as well as on archival IR data obtained with 
Spitzer and Herschel space observatories. The log of our ob¬ 
servations is presented in Table [U where the exposure time 
(Texp), the field of view (FOV), the final angular resolution 
(^), the final spectral resolution ((5A) and the bandwidth 
(FWHM) of the used filter are indicated. 
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Table 1. Summary of observational data 


Data set 

Filter or Disperser 

Date of obs 

Texp (sec) 

FOV 

en 

Sp. range 

SX or FWHM (A) 

Field #1 

FN655 

2014/10/16 

360 

6.1' X 6.1' 

3.5 

Hq+[Nii] 

97 


FN674 

2014/10/16 

360 

6.1' X 6.1' 

3.5 

[So] ^ 

60 


FN641 

2014/10/16 

160 

6.1' X 6.1' 

3.5 

continuum 

179 


FN712 

2014/10/16 

160 

6.1' X 6.1' 

3.5 

continuum 

209 

Field #2 

FN655 

2014/10/16 

180 

6.1' X 6.1' 

4.1 

Hq+[Nii] 

97 


FN674 

2014/10/16 

180 

6.1' X 6.1' 

4.1 

[So] ^ 

60 


FN641 

2014/10/16 

80 

6.1' X 6.1' 

4.1 

continuum 

179 


FN712 

2014/10/16 

80 

6.1' X 6.1' 

4.1 

continuum 

209 

long-slit #1 

VPHG1200@540 

2014/10/16 

1800 

6.1' X 1.0" 

2.6 

3690-7280A 

4.5 

long-slit #2 

VPHG1200G 

2014/08/03 

1800 

6.1' X 1.0" 

1.3 

3690-5680A 

5.4 

long-slit #3 

VPHG550G 

2014/08/22 

1800 

6.1' X 0.75" 

1.9 

3710-7890A 

9.5 

long-slit #4 

VPHG1200G 

2014/08/03 

1800 

6.1' X 1.0" 

1.4 

3690-5680A 

5.4 

long-slit #5 

VPHG1200G 

2014/07/31 

870 

6.1' X 1.0" 

1.5 

3690-5680A 

5.5 

FPI 


1990-2001 

1800 

10' X 10' 

3-4 

Hq 

0.2-0.35 


2.1 Spectral observations 


The spectral observations were performed with the 6-m tele¬ 
scope using the multi-mode SCORPIO (Spectral Camera 
with Optical Reducer for Photo metrical and Interferometri - 
cal Observations) focal reducer (lAfanasie v MoiseevI 2005ll 
and its new SCORPIO-2 version (|Afanasiev &: Moiseev! 

operating in the long-slit spectrograph mode. The 
scale along the slit was 0.36 arcsec per pixel. The CCD detec¬ 
tors EEV 42-40 in SCORPIO and E2V 42-90 (with the low 
sensitivity at the blue end of the spectrum) in SCORPIO- 
2 were employed. We observed the region around vdB 130 
with five slit positions. The information on all the obtained 
spectra is summarized in Table [H 

Data reduction was performed in a standard way us¬ 
ing the IDL software package developed at the SAO RAS for 
reducing long-slit spectroscopic data obtained with SCOR¬ 
PIO, which includes bias subtraction, line curvature and flat 
field corrections, linearization, and air-glow lines subtrac¬ 
tion. Dispersion correction of the spectra was performed us¬ 
ing the reference spectrum of a He-Ne-Ar lamp obtained 
during the observations. 

The observed region is densely filled with nebular emis¬ 
sion (see Section oil . Because of that we were forced to 
construct the night-sky spectrum model using only a few po¬ 
sitions along the slit in the direction of the dark cloud where 
the contribution of nebular emission to the spectrum is less 
significant. We interpolated these values of air-glow line in¬ 
tensities for all positions along the slit. The slit #3 does not 
cross the dark cloud, so in order to subtract the air-glow 
lines, we observed a nearby blank field where the nebular 
emission is absent. To convert the spectra to the absolute in¬ 
tensity scale, we observed the spectrophotometric standards 
BD+25d4655, BD+28d4211 and BD+33d2642 immediately 
after the object at a close zenith distance. Gaussian fitting 
was applied to measure the integrated fluxes of emission 
lines. Eor this purpose we used our own IDL software adopted 
for SCORPIO data. In order to estimate uncertainties of the 
gaussian approximation a number of synthetic spectra with 
predefined signal-to-noise ratios were generated and fitted 
with the IDL procedure. Then, the obtained uncertainty for 
a particular signal-to-noise ratio was used as an estimate of 
the uncertainty of the method. 

To estimate the Hll line-of-sight velocities in the re¬ 


gion, we used the results of long-term (1990-2001) Ha 
observations of the Cygnus gas-dust complex performed 
with the Eabry-Perot interferometer (EPI) mounted in the 
Cassegrain focus of t he 125-cm reflector of the Crim ean Lab¬ 
oratory of SAI MSU ([Lozinskaya et al1ll997l . [l998l L The pa¬ 
rameters of the instruments are listed in Table m We fitted 
the line profiles by one or several Gaussian curves assum¬ 
ing that the halfwidth of each profile component should be 
greater than that of the instrumental contour and signal- 
to-noise ratio is no less than 5. The peak velocities of each 
profile component were determined. 


2.2 Narrow-band images 

Optical images of two overlapping fields around vdB 130 
were taken at the primary focus of the 6-m telescope 
of SAO RAS with SCORP IO-2 multi-mode focal reducer 
([Afanasiev Moiseevll201ll i using filters EN655 and EN674, 
corresponding to Ha and [Sli] emission lines with central 
wavelengths of 6559 and 6733 A and EWHM = 97 and 60 A, 
respectively. The resolution was 0.36 arcsec per pixel. Since 
the EWHM of the used EN655 filter is broader than the sepa¬ 
ration between Ha and [Nil] emission lines, the image in this 
filter is contaminated by [Nil] 6548, 6584 A emission. Given 
the results of our spectral observations (see Section 14.21) . 
and taking into account the EN655 filter transmission, we 
conclude that the contribution of [Nil] lines to the obtained 
image may be up to 0.4 of the pure Ha emission. Hereafter 
we indicate with Ha the emission in the Ha line with this 
additional contribution due to the [Nil] lines. 

We used broader EN641 and EN712 filters centered on 
the continuum near the Ha and [Sil] emission lines to sub¬ 
tract the stellar contamination from the images obtained 
during the same night. In order to calibrate the emission 
line images to energy fluxes we observed the standard star 
BD+25d4655 immediately after observing the cluster region. 

After the data reduction has been made in IDL in the 
standard way (bias subtraction, flat field normalization, cos¬ 
mic ray rejection, background estimation and subtraction) 
we combined images of both observed fields around vdB 130 
in mosaics in [Sli] and Ha lines using our own iDL procedure 
and assuming positions of several stars in one frame as a 
reference. 
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2.3 Archival data 

This work is partially based on data from Spitzer and Her- 
schel space telescopes. The region studied was observed 
by Spitzer space telescope within the framework of the 
‘A Spi tzer Legacy Surve y of the Cygnus-X Complex’ pro¬ 
gram llHora et all 1201 ih . The images at 3.6/xm, 4.5/xm, 
5.8 /xm, 8.0 /xm, and 24 /xm were downloaded from the project 
web-pag43 at the NASA/IPAC Infrared Science Archive. 
The Herschel archival data obtained within the framework 
of the ‘Hi-GAL : The Herschel Infra red Galactic Plane Sur- 
veyO program (|Molinari et cTI[ 20101 1 were downloaded from 
the Herschel Science Archivqj. 

In thi s study we also used images from the 2MASS 
IR survey (|Skrutskie et al.l[2006l l and stellar proper mot ions 
data from the UCAC4 catalogue (jZacharias et al.ll2013l b 


3 OPEN STAR CLUSTER VDB 130 

Initially designation ‘vdB 130’ has been applied to a re¬ 
flection nebula associated with the BD+38 3993 star (HD 
228789, vdB 130a) of spectral class B1-B2 (|van den BerghI 
1 19661 : [Racind [l968l b The reflection nebula vdB 130 can be 
seen o n the optical maps shown in Figures [1^ and [Racind 
(|l974l l used photometric and spectroscopic data to select the 
group of 14 stars within 6^ from the nominal cluster centre, 
which could be considered as possible members of the sparse 
embedded cluster vdB 130 (Fig. (2). He noted an anomalous 
extinction law with AyIE{B — V) = 8.1 ±1.2, probably due 
to inhomogeneity of the interstellar medium in the proxim¬ 
ity of the cluster. I Kharchenko et ahl (|2013l L with pipe-line 
processing technique, estimate d the cluster parame ters using 
data from PPMXL catalogue (jRoeser et al.ll20iq j and con¬ 
sidering stars with proper motion components close to the 
centroid with (/Xa,/X 5 ) ~ (—2.22,—4.39) masyr“^ as cluster 
members. From 2MASS JHKs photometric data they de¬ 
rived a cluster age log(t) ~ 7.4, distance d ~ 1.6 kpc, colour 
excess E(J — H) ^ 0.29 mag, infrared interstellar absorption 
Aks ~ 0.30 mag, and apparent radius of about 6b 


3.1 Revision of the stellar content of vdB 130 

Most possible members of vdB 130 move parallel to 
the galactic plane wi th mean proper inotion s {iii^fih) ~ 
(—4.9, —0.6) masyr“^ ([Kharchenko et al]l2013l b In order to 
estimate the contribution of the disk differential rotation 
and that of the Sun peculiar velocity to the observed proper 
motions of cluster stars for heliocentric distances ruei ~ 
1.5 — 2.5 kpc, the parameters of Milky Way disk rotation 
and Sun’s velocity components with respect to the LSR 
(Local Standard of Rest) were taken from a detailed kine- 
matical analysis of young stell ar populations performed by 
IZabolotskikh et ^ (|2002l j and I Melnik fc Damibisl (|2009l L If 
we use a currently adopted value of the Sun’s distance 
from the Milky Way centre, 8.3 kpc, all helds located 
along the line of sight in Gygnus directions with helio¬ 
centric distances in the ruei ~ 1.5 — 2.5 kpc range, are 


^ https://www.cfa.harvard.edu/cygnusX/ 

^ https://hi-gal.iaps.inaf.it/higal/ 

^ http://www.cosmos.esa.int/web/herschel/science-archive 
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Figur e 2. Star cluster vdB 130. Stars from the original iRacinel 
(I 1974 I) paper are marked with indices from Ir to 14r. Stars se¬ 
lected in this study based on their pro per mot i ons a re marked 
with numbers from 1 to 36. Stars from iRacin i iIi974i list with 
proper motions outside of the adopted intervals are shown with 
red colour. The letter a marks the star K4 projected onto an 
extended IR source (Blob E; see section 4.3). The central star 
BD-h38 3993 of the reflection nebula is marked as llr. 


nearly equally distant from the galactic centre, ^ S.l kpc. 
From the 3D Bottlinger equations we estimate mean appar¬ 
ent proper motion at 1.5-1.8 kpc heliocentric distances as 
(/xz(rot), /j.b{rot)) ^ (—5.4,—0.9) masyr“^. 

These values are in a very good agreement with 
the proper motion com ponents for vdB 130 derived by 
I Kharchenko et all (|2013l L Therefore, we c annot exclude that 
stars assigned to the open cluster by iKharchenko et ^ 
(|2013h represent in fact a non-bound stellar population of 
the extended association Gyg OBI, participating in over¬ 
all galactic differential rotation. Moreover, the accuracy of 
available proper motions in this held does not allow distin¬ 
guishing reliably true members of an unbound OB associa¬ 
tion from_Joregroimb_anb_b^kgromd held stars. 

iLe Duigou KnodlsedeJ ( 2002h also used 2MASS data 
to study the stellar population and the structure of vdB 130, 
along with other 21 clusters in the Gygnus region. Un¬ 
fortunately, the lack of detailed data on this cluster does 
not allow estimating its parameters reliably. It shou ld be 
noted that application of King’s empirical density law ([Kind 
Il962h to describe the apparent distribution of stars in 
dense and highly non-uniform stellar helds distorted by in¬ 
homogeneous absorption is unlikely able to provide reli- 
able estimates of import a nt str ucture cluster parameters. 
iLe Duigou Knodlsederl (|2002h used the cluster distance 
derived bv [Racing (1973) and claimed vdB 130 as the richest 
cluster among those studied by them in the Gygnus region. 
It should be noted, however, that they used no kinematical 
criteria to select probable cluster members. 

Despite the insufficient evidence supporting the nature 
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Table 2. Candidate members of vdB 130 


Star # 

UCAC4 



BS 

Tycho-2 

2Mkey 

J (mag) 

H (ma 

g) 


K (mag) 


B (mag) 

V (mag) 

1 

647-086881 

-1.5 ± 6.0 

-5.7 

± 

5.9 


305146782 

13.2 


0.03 

12.856 

± 

0.03 

12.666 ± 

0.03 

16.031 

15.332 

2 

647-086916 

-5.4 ± 2.2 

-4.2 


2.2 


305146998 

12.01 


0.02 

11.679 

± 

0.02 

11.56 ± 

0.02 

14.212 

13.582 

2r 

647-086761 

-5.0 zb 2.2 

-2.5 


4.5 


305106035 

11.703 

± 

0.02 

11.352 


0.02 

11.225 ± 

0.02 

13.907 

13.28 

3 

647-086873 

-4.5 ± 1.9 

-1.6 

± 

2.4 


305105199 

12.045 


0.02 

11.95 


0.02 

11.853 ± 

0.02 

13.555 

13.067 

3r 

647-086871 

-4.8 zb 0.9 

-4.5 


0.6 

3151-00092-1 

305105710 

9.591 

zb 

0.02 

9.421 

zb 

0.02 

9.33 zb 

0.02 

11.597 

10.99 

4 

647-086802 

-3.5 zb 6.0 

-5.1 

± 

2.8 


305105443 

12.615 

± 

0.02 

12.33 

zb 

0.02 

12.203 ± 

0.02 

15.636 

14.883 

4r 

647-086813 

-4.4 zb 1.9 

-5.3 


3.7 


305105485 

11.376 

± 

0.02 

11.086 

± 

0.02 

10.909 ± 

0.02 

13.852 

13.125 

5 

647-086786 

-0.2 zb 4.5 

-3.9 

± 

4.3 


305105367 

13.047 

± 

0.02 

12.316 

± 

0.02 

12.164 ± 

0.02 



6 

647-086831 

-3.2 zb 4.8 

-5.8 

± 

4.5 


305105399 

12.948 

± 

0.02 

12.618 


0.02 

12.469 ± 

0.03 

15.634 

14.914 

7 

647-086791 

-2.9 zb 5.4 

-3.1 

± 

4.8 


305106069 

13.761 

± 

0.03 

13.408 

± 

0.02 

13.282 ± 

0.04 

16.781 

15.847 

8 

648-090623 

-3.1 zb 2.1 

-1.1 

± 

1.7 


305106265 

11.65 

zb 

0.04 

11.382 

± 

0.02 

11.32 zb 

0.02 

13.295 

12.727 

8r 

647-086860 

-6.1 ± 4.2 

-4.3 

± 

4.2 


305105753 

11.994 

± 

0.02 

11.4 


0.03 

10.962 ± 

0.02 

15.466 

14.642 

9 

648-090627 

-0.6 zb 1.6 

-5.2 

± 

2.7 


305106343 

11.635 

± 

0.02 

11.311 

± 

0.02 

11.239 ± 

0.02 

13.801 

13.09 

10 

648-090649 

-4.0 ± 2.2 

-5.9 

± 

2.0 


305106565 

12.028 


0.02 

11.75 


0.02 

11.676 ± 

0.02 

13.986 

13.352 

lOr 

647-086781 

-4.0 ± 3.3 

-3.9 

± 

4.0 


305105837 

11.176 

± 

0.02 

10.808 


0.02 

10.652 ± 

0.02 



11 

648-090655 

-3.8 ± 4.6 

-6.8 


4.4 


305106417 

13.168 

± 

0.02 

12.792 


0.02 

12.572 ± 

0.03 

16.971 

15.835 

Hr 

647-086776 

-4.3 ± 0.8 

-4.9 

± 

0.5 

3151-00165-1 

305105806 

8.98 


0.03 

8.816 


0.02 

8.745 ± 

0.02 

10.755 

10.214 

12 

648-090662 

-0.4 ± 5.9 

-5.3 


6.0 


305145860 

12.141 

± 

0.02 

11.543 

± 

0.02 

11.354 ± 

0.02 

16.503 

15.197 

12r 

647-086785 

-4.4 zb 0.9 

-5.4 


1.1 

3151-00109-1 

305105748 

9.79 

zb 

0.02 

9.621 

zb 

0.02 

9.526 zb 

0.02 

11.54 

11.022 

13 

648-090657 

-6.8 ± 5.2 

-3.6 


5.8 


305106346 

14.05 


0.02 

13.73 


0.04 

13.534 ± 

0.05 



13r 

647-086775 

-3.9 zb 4.4 

-7.0 


7.0 


305105714 

11.738 

± 

0.02 

11.379 

± 

0.02 

11.236 ± 

0.02 



14 

647-086865 

-6.1 zb 4.4 

-2.0 

± 

4.5 


305106215 

13.2 

zb 

0.02 

12.664 

± 

0.02 

12.457 ± 

0.03 

17.236 

15.873 

15 

647-086859 

-1.3 zb 5.5 

-6.4 


5.3 


305106104 

14.472 

± 

0.04 

13.829 

± 

0.04 

13.734 ± 

0.05 



16 

647-086890 

-2.1 zb 2.7 

-5.5 


3.1 


305146111 

13.146 

± 

0.03 

12.919 

± 

0.03 

12.786 ± 

0.03 

15.409 

14.717 

17 

647-086910 

-2.6 zb 2.3 

-4.5 

zb 

3.7 


305146062 

12.202 

± 

0.02 

11.907 

± 

0.02 

11.829 ± 

0.02 

14.278 

13.713 

18 

647-086874 

-1.4 ± 2.6 

-3.8 


2.4 


305105912 

11.632 

± 

0.02 

11.385 


0.02 

11.233 ± 

0.02 

13.898 

13.214 

19 

647-086887 

-6.3 zb 0.9 

-3.6 

± 

1.1 


305146507 

11.091 

± 

0.02 

10.943 

± 

0.02 

10.839 ± 

0.02 

13.003 

12.406 

20 

647-086833 

-0.2 ± 4.6 

-4.6 

± 

4.4 


305105567 

13.154 

± 

0.03 

12.775 


0.03 

12.58 ± 

0.03 

16.57 

15.55 

21 

647-086826 

0.4 ± 3.5 

-5.5 


2.4 


305105664 

12.03 

zb 

0.03 

11.581 

± 

0.04 

11.402 ± 

0.03 

15.08 

14.146 

22 

647-086829 

-0.8 ± 4.3 

-6.5 


4.3 


305105800 

12.741 

± 

0.99 

12.166 

± 

0.03 

11.925 ± 

0.02 

16.35 

15.262 

23 

647-086940 

-2.4 zb 6.3 

-5.7 


6.1 


305146927 

12.128 

± 

0.02 

11.252 

± 

0.02 

10.956 ± 

0.02 



24 

647-086954 

-5.6 zb 4.9 

-2.4 


1.8 


305146496 

12.121 

± 

0.02 

11.708 

± 

0.02 

11.457 ± 

0.02 

14.536 

13.87 

25 

647-086949 

-3.1 zb 6.5 

-6.5 


6.2 


305146341 

12.857 

± 

0.03 

12.081 

± 

0.03 

11.835 ± 

0.03 



26 

648-090681 

-0.9 zb 7.2 

-2.2 


6.8 


305145870 

12.397 

± 

0.02 

11.467 

zb 

0.02 

11.183 ± 

0.02 



27 

648-090678 

0.2 ± 3.9 

-4.5 


4.8 


305145626 

13.078 

± 

0.02 

12.736 

± 

0.02 

12.655 ± 

0.03 

15.483 

14.753 

28 

648-090639 

-5.3 ± 2.2 

-3.4 


2.4 


305106466 

12.516 

± 

0.02 

12.141 

± 

0.02 

12.072 ± 

0.02 

14.73 

13.959 

29 

648-090628 

-3.1 zb 4.7 

-5.6 

± 

4.4 


305106531 

13.004 

± 

0.02 

12.611 

± 

0.02 

12.463 ± 

0.02 

15.776 

14.911 

30 

647-086696 

-1.1 ± 4.5 

-2.6 

± 

4.4 


305036100 

13.573 

± 

0.02 

13.161 

± 

0.02 

13.095 ± 

0.03 

16.62 

15.574 

31 

647-086727 

-0.1 ± 4.6 

-3.0 

± 

4.6 


305105507 

13.653 


0.03 

13.263 

± 

0.04 

13.203 ± 

0.04 

16.397 

15.427 

32 

647-086773 

-5.9 ± 4.3 

-5.8 

± 

4.3 


305105002 

12.236 

± 

0.02 

11.49 


0.02 

11.275 ± 

0.02 



33 

647-086823 

-1.4 ± 4.8 

-7.0 

± 

4.6 


305104995 

12.421 

± 

0.02 

11.567 

± 

0.02 

11.295 ± 

0.02 



34 

647-086856 

-3.5 ± 4.0 

-7.0 

± 

1.3 


305105078 

12.536 

± 

0.02 

12.221 

± 

0.02 

12.074 ± 

0.02 

15.613 

14.701 

35 

647-086866 

-0.2 ± 5.9 

-2.9 

± 

6.3 


305105086 

13.673 

± 

0.03 

13.354 

± 

0.03 

13.144 ± 

0.04 

17.342 

15.927 

36 

647-086897 

-1.7 ± 5.8 

-7.8 


6.1 


305146976 

13.636 

± 

0.03 

13.261 

± 

0.03 

13.086 ± 

0.04 

16.644 

15.926 

Ir 

647-086816 

23.9 ± 18.8 

-27.0 ± 9.0 


305105825 

11.594 ± 0.03 

10.744 ± 0.02 

10.401 ± 0.02 

15.085 

14.225 

5r 

647-086841 

-7.0 ± 4.4 

2.9 ± 4.4 


305105774 

11.988 ± 0.04 

11.435 ± 0.04 

11.108 ± 0.04 

15.744 

14.741 

6r 

647-086746 

53.1 ± 2.8 

34.5 ± 1.8 


305105807 

10.211 ± 0.02 

9.666 ± 0.02 

9.56 ± 0.02 

13.152 

12.124 

7r 

647-086750 

32.4 ± 5.7 

-8.5 ± 6.1 


305105772 

13.432 ± 0.02 

12.924 ± 0.03 

12.835 ± 0.03 



9r 

647-086849 

-4.7 ± 1.7 

-9.0 ± 3.3 


305105989 

11.330 ± 0.02 

11.092 ± 0.02 

11.025 ± 0.02 

13.383 

12.752 

14r 

647-086779 






305105685 

12.536 ± 0.02 

12.163 ± 0.02 

11.945 ± 0.02 




of vdB 130 as a true gravitationally bound stellar cluster, 
further study of the stellar content of this star field, located 
at the distance of about 30-40 parsecs from the centre of 
the Cyg OBI association, seems to be worthwhile. To find 
new possible members of vdB 130 and confirm th e mem¬ 
bershi p of stars earlier attributed to this cluster by iRacind 
illaii) we reconsidered the data of 175 stars in a rather ex¬ 
tended region with the diameter of 12 arcmin and central co¬ 
ordinates (a 2 ooo ~ 20^ 17”^50®, fcooo ~ 39°21'25", for which 
stellar proper motions can be taken from UCAC4 catalogue 
(|Zacharias et al.l l2013l ) . The concentration of st ars around 
the centroid dehned bv I Kharchenko et al.l (|2013l i is evident 
in the top panel of Fig. O This can be readily explained by 
the dominant contribution of disk differential rotation and 
solar motion to the apex. 

We selected 36 UCAC4 stars as possible members of the 
cluster, with proper motion vectors that differ from the con¬ 
sidered centroid by less than 4masyr“^ (see bottom panel 
of Fig. [3}. This choice of a limiting value corresponds to 


the typical proper motion error in the UCAC4 catalogue, 
while in most cases the errors range from 1 to 10 mas yr“^. 
Adopting a larger threshold of proper motions, we risk of 
contaminating our sample with a large number of foreground 
and background field stars with large transversal velocities. 
(Note that the difference of 4masyr“^ at the distance of 
1.5 kpc is equivalent to 28kms“^ difference in transversal 
velocity. This value greatly exceeds the typical velocity dis¬ 
persion in OB associations.) 

Eight of 14 stars from iRacin i (Il974h original list sat¬ 
isfy our proper motions criteria (Fig. [2]). The extended list 
of possible cluster members is given in Table [2l Star iden¬ 
tifications, proper motions and respective errors, 2MASS 
identifications and JHKs magnitudes, BV magnitudes (for 
some stars) were taken from the UCAC4 catalogue and are 
printed keeping original catalogue format. T he colu i nns o f 
Table [2] contain the star number according to iRaciii idliTi) 
list (with a letter ‘r’) and/or our list, its identification in the 
UCAC4 catalogue, proper motion components /Xa, iJis and 
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fia, mas/yr 


Figure 3. (Top) Proper motion (iia, f^s) diagram for stars from 
the UCAC4 catalogue in the circle within 6' from the centre of 
the field. (Bottom) Central part of the diagram. Filled circles 
indicate possible cluster members selected in this study. Shown 
with triangles are stars with the most accurate proper motions 
taken from the TYCH02 catalogue. Field stars are shown with 
open circles. 


corresponding errors, TYCHO2 catalogue (|Hog et al.ll2002l ) 
identification (for some stars), 2MASS catalogue identifica¬ 
tion (2Mkey) and JHKs magnitudes and errors, and BV 
magnitudes. In the bottom part of the table we show the 
data for stars from iRacin 3 (Il974l) list with motions outside 
of the adopt ed inter v als (I r, 5r, 6r, 7r, 9r and 14r). Stars Ir 
and 5r from [Racinel (| 19741 ) original list, which are formally 
not considered as the cluster members, apparently have large 
velocities and move toward each other, with velocity compo¬ 
nents relative to vdB 130 (V^, 14) ~ (—31 zb 137, —260 zb 83) 
and (26 zb 47,61 zb 9) kms“^, respectively. However, their 
proper motions are measured with very large errors, partly 
due to their asymmetric shapes in the optical images. There¬ 
fore, we cannot decisively exclude these two bright stars from 
the vdB 130 cluster members. 

Based on the sample of stars shown in Table [2l we 
estimate basic characteristics (colour excess and distance) 
vdB 130, using JHKs and BV magnitudes and isochrone 
fitting technique. A colour-absolute magnitude diagram 




Figure 4. Colour-absolute magnitude diagrams My vs. {B — V) 
(top) and Mks vs. ( J — H) (bottom) for vdB 130 possible mem¬ 
bers selected by proper motions (circles). Theoretical isochrones 
with ages t = 5,10, 30,1 00, 300 Myrs (frorn top to bottom) are 
also s hown for BV data (iGirardi et ^ ]l2nn2h and 2MASS JHK. 
data llBonatto et al.ll20n4h . Diamonds indicate stars from lRacinj 
lIlQTdh list with large differences in proper motions from the cen¬ 
troid. 


(CMD) V vs. B — V for 34 probable vdB 130 members with 
avail able BV data and Padova isochrones (|Girardi et al.l 
I2OO2I) calculated for solar abundance and for ages t = 
5,10, 30,100, 300 Myrs are shown on the top panel of Fig. [H 
If we take into account typical photometric uncertainty of 
0.03 — 0.05 m ag and the real overdens ity of bright stars, 
clearly seen on iKharchenko et al.l (|2013l ) atlas for vdB 130, 
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Figure 5. A (/ia, diagram for stars with most precise mea¬ 
surements of proper motions. 


the eye-fit gives an estimate for the cluster age of about 5- 
10 Myrs, assuming that the best-fit isochrone is a lower left 
envelope for stars with lowest colour indices. 

Despite a relatively small scatter in proper motions, 
colour scatter along the {B — V) axis is quite large (up to 
0.6 mag for faint stars), confirming the substantial differen¬ 
tial reddening over the cluster area. Minimal colour excess 
E{B — V) = 0.79 zb 0.02 mag, and the apparent distance 
modulus {V — Mv)app = 15.0 zb 0.3 mag. Under the as¬ 
sumption of a ‘normal’ extinction law (|Cardelli et al.|[l989l : 
iMartin Whitt^ll990l b with Rv = AyIE(B — y) 3.08, 
true distance modulus can be estimated as (y — Mv)o ^ 
12.6 zb 0.3 mag. In this case the distance to vdB 130 
stars having minimal redden ing values can be e s timat ed as 
3.3 zb 0 . 4 kpc , in excess of iKharchenko et all (|2013l i and 
iRacind (|l974l i estimates. Of course, our estimates of the 
colour excess and the distance relate to bluest and brightest 
stars on the CMDs, which are expected to li e on the front 
side o f vdB 130. It should also be noted that ISchlegel et al.l 
(|l998h and ISchlaflv Finkbeinerl (|201lh give total line-of- 
sight colour excess E{B — V) values in ranges 4.85 — 9.6 mag 
and 4.2 — 8.3 mag, respectively, for a 5 arcmin radius cir¬ 
cle, centred in the field of the vdB 130 cluster. These data 
support Racine’s conclusion on possible large differential ex¬ 
tinction in this region. 

A somewhat smaller distance can be derived from in¬ 
frared data. A Ks vs. (J — H) colo ur-magnitude diagra m for 
36 stars and the isochrones from (|Bonatto et al.|[2003 b cal¬ 
culated for the same ages as for optics, are shown on the bot¬ 
tom panel of Fig.|4l Using the same fitting technique, we esti¬ 
mate colour excess E{J—H) ^ 0.27zb0.02 mag and apparent 
distance modulus {K — MKs)app ~ 11.57 zb 0.3 mag. Under 
the assumption of a ‘normal’ extinction law (|Dutra et al.l 
I2OO3I : iHe et al]ll995l b with Aks ~ 1.13E'(J — if), we esti¬ 


mate near infrared (NIR) absorption as Aks ~ 0.31 zb 0.03, 
true distance modulus as {K—Mks)q ~ 11.26zb0.3 mag, and 
the corresponding cluster distanc e as 1.8 zb 0.3 kpc, in a gen - 
eral agreement with the result of IKharchenko et al.l (|2013l b 

The differences between true distance moduli estimated 
independently from optical and NIR data can be attributed 
to abnormal features of the extinction law in the cluster 
direction. Obviously, any deviations from the generally ac¬ 
cepted ‘normal’ extinction law should manifest themselves 
much stronger in the optical range than in the infrared. To 
match the distance estimates obtained from optical and in¬ 
frared data, i.e. from V vs. (B — V) and Kg vs. (J — if) 
CMDs, we need to assume that Ry in this field differs 
from the ‘standard’ va lue of 3.08 and reaches 4 .5-4.8 , in 
a good agreement with iMoffat Schmidt-Kalerl (|l976l i re¬ 
sult. Note, however, that this conclusion is rather prelimi¬ 
nary, because distance moduli in optical and infrared have 
been really estimated from only 6-7 stars with the lowest 
colour excess. Moreover, some of these stars are located in 
the upper left part of the CMD, where the slope of the main 
sequence and isochrones is large enough and, thus, vertical 
shift of isochrones (i.e. apparent distance modulus) is deter¬ 
mined not as reliably as the colour excess. 

We now consider the group of stars with the most ac¬ 
curate me asurements of pr oper motions. Three stars from 
TYCH02 (|Hog et al.ll2002h are notable among them, with 
proper motion errors smaller than l.lmasyr“^. On the 
proper motion diagram seven stars form a compact 

group with 4masyr“^ diameter (Fig. [5|). Four of them, 3r 
(TYC 3151-00092-1), Hr (TYC 3151-00165-1), 12r (TYC 
3151-00109-1), and 19 (2MASS 305146507), are included 
in our list and they are the brightest stars in our sample. 
Mean components of their proper motion are approximately 
= (—4.95, —4.60) masyr“^ with the RMS scatter of 
(0.9, 0.8 ) masyr~^. The differenc e between the mean motion 
and the IKharchenko et al.l (|2013h centroid (and, at the same 
time, contribution of disk differential rotation) is approxi¬ 
mately 2.7zb0.9 masyr“^. At the cluster distance (~ 1.8 kpc) 
this leads to 20kms“^ difference of transversal velocities. 
We can conclude that these brightest members of the young 
open cluster vdB 130 associated with the Cyg OBI give an 
idea on the apparent cluster motion on the sky. To draw 
more definite conclusions on the reality of the gravitation¬ 
ally bound young stellar group additional precise proper 
motions measurements and spectroscopic observations are 
needed that would allow refining the spectral classification 
and provide accurate radial velocity measurements. 

3.2 Spectra of vdB 130 stars 

We used the spectra described in Section 2.1 to determine 
the spectral types of the stars Ir, 4r, 5r, 6r, Hr, and 22 
(Table [2} and of the a star (Fig. [2}. All the stars, except 
for the last one, have been classified as B stars. We anal¬ 
ysed the a star in order to assess its possible contribution 
to the emission of the gas-dust nebula studied in Section 4.3 
and classified it as a K4 star. Note that the spectral energy 
distribution of this star is practically undistorted by inter¬ 
stellar reddening, which is unlikely for a star located at the 
distance of vdB 130 (1.5-2 kpc), and this further confirms 
that a is a foreground field star. 

According to our estimates, the spectral type of Ir and 
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Hr is BIV up to 0.5 subtype; that of 4r and 5r is interme¬ 
diate between BIV and B2V; and that of 6r and 22 is B5- 
B6V with an accuracy of 1 subtype. (We adopted the spec¬ 
tra of MK standards from I Gray CorballvlEoOQl .) All the 
obtain ed spectra are shown in Fig. [G] According to iRacind 
(|l974l b the spectral types of stars Ir, 5r, and Hr are B1-2V, 
B3-5V and B2III, respectively. However, in spectra of B3- 
5V-type sta rs conspicuous MpII 4 481 A absorption should 
be present (|Grav V Corball^ l2009l b while it is barely seen 
in the spectrum of star 5r (Fig.[6|) that was obtained with a 
substantiall y higher disper sion than the spectrum obtained 
byE acin i (HOT i)M acinel (|l974l l estimated the luminosity 
class of the star Hr as III based on smaller widths of its 
Balmer absorption lines. However, the widths of these lines 
in our spectra of stars Hr and 5r are nearly the same (Fig. 
El). Furthermore, one should take into account the presence 
of Balmer emission lines in the gas-dust nebulae surrounding 
the considered stars which can distort the actual appearance 
of stellar absorption lines. 

Approximating the total interstellar and circumstellar 
extinction by the normal interstellar extinction law allows 
estimating the colour excess E(B — V) from the spectral 
energy distribution (SED). The colour excesses of the inves¬ 
tigated stars differ slightly from each other, being equal to 
1.1 mag for 4r and Hr, 1.3 mag for Ir, 1.4 mag for 5r and 
22, and 1.5 mag for 6r. The relative continuum flux differ¬ 
ence between the dereddened SEDs of all the studied stars 
and th e SEDs of standard stars of the correspon ding spectral 
types (jValdes et al.ll200j : ISilva V Cornelll[r992l l does not ex¬ 
ceed 10 per cent at all wavelengths in the considered spectral 
band. 

Red spectra (up to ~ 7300 A) were taken only for stars 
Ir and 4r. In the case of the star Ir the red spectrum com¬ 
plements rather well the spectrum covering the 4400-5580 
A wavelength interval and yields the extinction parameters 
which do not differ substantially from those inferred from 
the bluer spectrum. On the other hand, the continuum of the 
4r star is reproduced much better assuming a non-standard 
total-to-selective extinction ratio of i^v = 4.0, which agrees 
with the estimate presented in the previous section. The 
colour excess in this case \s E{B — V) — 1.05 mag. 


4 INTERSTELLAR MEDIUM IN THE VDB 130 
AREA 

4.1 Molecular cloud and ionized hydrogen in the 
vicinity of the cluster 

VdB 130 is located in projection near the north¬ 
western border of an extende d optical, radio and IR 
shell around Gyg OBI (Fig. [H (iLozinskava V Sitnil3 Il988l : 
iLozinskava V RepinI Il990l : Saken et al.l Il992l b This super¬ 
shell is also clearly discernible at all wavebands o bserved 
with the Spitzer space tele scope (iHora et al ]l200flh and at 
large-scale CO distribution (ISchneid er et al.ll 2007lb The dis¬ 
tance to the association is 1.5 kpc (|Garmanv V Stencelll 19921 : 
ISitnik V Merniklll996h . From the analysis of the cluster stel¬ 
lar population, we estimate that within errors Gyg OBI and 
vdB 130 are located at the same distance (see sect. 3.1). 

The studied region is projected onto a comet ary molec¬ 
ular cloud (Fig. [T]3, [7^) having a length of ~ 0.4° (cloud 



Wavelength (A) 


Figure 6. Spectra of some vdB 130 stars corrected for the in¬ 
terstellar reddening and a standard B1.5V spectrum. Star desig¬ 
nations are the same as in Fig. [2] The inset shows an expanded 
view of a small part of the star 5r spectrum and the standard 
B3V spectrum. Main absorption lines are labelled. Shaded verti¬ 
cal bands show wavebands with incorrect background subtraction. 
Spectra are given in arbitrary units as we do not use an absolute 
calibration. 

A in I Schneider et al.l l2007l j . The cloud is stretched along 
the galactic latitude, with its ‘head’ protruding inside the 
Gyg OBI supershell. The shape of the GO cloud hints that 
its ‘head’ is actually a pillar on the inner supershell bor¬ 
der. Some stars of vdB 130 are seen in the region of the 
cloud head, and all the other stars surround it (Figs. [2l [7^). 
The GO emission from the cloud spans the v elocity interval 
Hlsr ~ — l-r5kms“^ ([Schneider et ahlboOTl b The fact that 
the median line-of sight velocity for 34 out of 70 Gyg OBI 
stars is 4km s~^ with a standard deviation of 8.9kms“^ 
([Sitnikl [2003l j suggests that the association and the cloud 
are close neighbours in space. The presence of a reffection 
nebula and considerable scatter of interstellar extinction es¬ 
timates toward cluster stars can also be an ind irect indicato r 
of a link between vdB 130 and the CO cloud (|Racinelll974l b 
The differential reddening could also be possible if the 
cluster is in the background of the cloud, but the extinction 
values found above are too low for this to be true. Indeed 
the minimum E(B — V) in the cluster area is 0.79 mag, 
which corresponds to Ay only of the order of 2.5 mag. On 
the other hand, some stars, projecting onto the cloud head, 
have colour excesses of the order of 1.3-1.5 mag. Assuming 
Ry ^ 4 _ 5 (see Section 3), the extinction to ward these stars 

S blly agrees with the estimates from I Schneider et al.l 
for their Cloud A (5-8 mag). 

Our FPI observations can be used to find kinematic 
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Figure 7. VdB 130 region with contours of the ‘head’ of the CO cloud (|Schneider et al.lboOTh superimposed onto the E-DSS2 map (a) 
and onto the 8.0/rm image (b). 


evidence in favour of the physical connection between the 
molecular cloud and ionized hydrogen in the cluster out- 
skirt. The data exhibit single-component Ha profiles toward 
vdB 130 with peak line-of-sight velocities ranging from 0 
to 15kms“^. These velocities characterize gas whose mo¬ 
tions are mostly due to Galactic rotation and streaming 
motions caused by spiral density waves. To analyse char¬ 
acteristic systematic motions of ionized hydrogen in the 
studied region, we use a ‘ diagram of occurrence ’ of line- 
of-sight velocities (see, e.g. I Arkhip ova et al.l l2013l and ref¬ 
erences therein). Characteristic line-of-sight velocities of the 
ionized gas Vlsr in the region are defined as average val¬ 
ues at half-maximum of the corresponding diagram together 
with root-mean-square errors of their estimates. The ob¬ 
tained interval of the most commonly occurring ionized gas 
velocities toward vdB 130 is AVlsr ~ 3 — 9kms“^ (based 
on 24 velocity measurements out of 26) and is shifted in 
the positive direction relative to velocities of CO emission 
Vl sk ^ — 1 T 2.5kms~^ in the head of the cloud (see Fig. 2 
in ISchneider et aI1l2007l i. 

Note that in the tail of the comet ary molecular cloud 
(about 5—12 arcmin westward of vdB 130), where the influ¬ 
ence of cluster stars should be less pronounced, Hii velocities 
(Hlsr ~ 3 —6kms“^) do indeed coincide with CO velocities 
in the same area (Vlsr ~ 2 — 5kms“^). This is an evidence 
in favour of the physical connection between the CO cloud 
and ionized hydrogen in its environment. Minor differences 
in Vlsr of molecular and ionized gas, observed in the cloud 
head and in its vicinity, may indicate the influence of winds 
and radiation from vdB 130 stars onto the surrounding gas. 

Inside the supershell around Cyg OBI, profiles with 
multiple peaks are observed. Specifically, in addition to the 
main component, they exhibit weak high-velocity features 
in the wings of the Ha line. These shifted components re¬ 
veal local high-velocity gas motions caused by stellar winds 
and/or supernova explosions (with a possible contribution 



Figure 8. Position-velocity diagram for the area shown by the 
white rectangle in Fig.[T] The star HD 193595 is located at the co¬ 
ordinate origin, and the molecular cloud is located at the distance 
of 0.°4 from the star. 


from distant Hll regions in the Perseus arm at large nega¬ 
tive velocities). 

The molecular cloud can be affected by UV radia¬ 
tion and stellar winds from at least two association mem¬ 
bers, namely, HIP100173 (07e) and HD193595 (08), lo¬ 
cated approximately along the extension of the cometary 
cloud axis, in 16 arcmin (7 pc) and 27 arcmin (12 pc) from 
the considered region (Fig. [T]). These stars are the near¬ 
est membe rs of the Cyg OBI associat i on with strong stel¬ 
lar w inds (jBlaha Humphrevsl Il989l : ICarmanv S tencell 
Il992l i. However, the entire vdB 130 star-gas-dust complex 
is located close to the boundary of the supershell surround¬ 
ing Cyg OBI. Therefore, the projected expansion velocity of 
the shell should be close to zero in the considered region. In- 
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Figure 9. Maps of ionized gas in the direction of vdB 130. Bright stars are masked in panels a) and c). Numbers indicate vdB 130 stars, 
(a) An Ha image and slit positions used for our spectral observations, (b) A false-colour RGB-image where red colour corresponds to 
[Sll] emission, green colour corresponds to Ha emission, and blue colour indicates the flux distribution in the Spitzer 8 /rm band, (c) A 
map of [Sll]6717,6731/Ha flux ratio with overlaid isophotes in S/rm band, (d) A map of uncertainties of the [Sll]6717,6731/Ha flux ratio 
shown in the panel (c) with overlaid isophotes in 8 p,m band. 


deed, as we pointed out above, spectra taken in the vicinity 
of vdB 130 exhibit mostly single-peaked Ha profiles. 

In Fig. [8] we show the PV-diagram, i.e., the plot of 
Hli line-of-sight velocities vs. the projected distance to the 
HD 193595 star. (The area for which this PV-diagram is con¬ 
structed is shown with a white rectangle in Fig. [1^.) High- 
velocity Hll motions are evident within 0.°4 from HD193595. 
This separation corresponds to the distance from the star to 
the outer boundary of the supershell and to vdB 130. 

Thus, the cometary shape of the cloud, the close spatial 
location of the cloud and the cluster, highly non-uniform 
reddening of vdB 130 stars, and the distribution of line-of- 
sight velocities of ionized gas in the cloud neighbourhood 
are all indicative of a possible physical connection between 
the CO cloud, Cyg OBI, and vdB 130. 

However, the exact nature of this connection is not im¬ 
mediately clear. CO emission implies relatively high density 
and correspondingly high extinction, so that some cluster 
stars are at least partially embedded in the absorbing ma¬ 


terial with Ay of the order of 3-6 mag. However, the CO 
emission in tensity in the directio n of vdB 130 implies higher 
extinction ([Schneider et al.|[2007l i (see above) so that at least 
some CO-emitting gas should be located behind the cluster. 
It is also possible that some CO emission comes from the 
dense material related to infrared features discussed below. 


4.2 Emission spectrum of the interstellar medium 
toward vdB 130 

Figure [9] shows various [Sli] 6717,6731 and Ha emission maps 
of the area surrounding vdB 130, based on our observations 
performed with the 6-m telescope of the SAO RAS. An Ha 
image is shown in Fig. [9^. In Fig.fSJ) the Ha map is combined 
with a [Sll] map and an 8 /im Spitzer map. 

In the interstellar medium surrounding stars lOr, Hr, 
12r, and 13r a well-defined stratification of Ha and [Sil] emis¬ 
sion can be seen, which is typical of Hii regions. Specifically, 
a ‘green region’ dominated by Ha emission surrounds the 
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Figure 10. Profiles of the [Oiii]5007/H/3, [Nii]6548,6584/Ha, and [Sii]6717,6731/Hq; flux ratios along four slits shown in Fig. [9^-. The 
filled grey areas indicate the uncertainties of the measured flux ratios. The [Olll]5007/H/3 ratio is shown by the lighter shade of gray 
and that of the [Sii]6717,6731/Hq; ratio is shown by the darker shade. The larger filled circles indicate the regions where we infer the 
combined emission mechanisms (see text). 


stars mentioned above in the centre and the lower part of 
Fig.[9l3. (Blue halos around stars Fig. [91: are wings of a point 
spread function in Ha). In the north-western part a ‘yellow 
region’ is located, where both Ha and [Sll] lines are bright. 
(Note that a similar region in the south may be partly ob¬ 
scured by the material responsible for the cluster reddening). 
The periphery of the Ha emitting area is dominated by ‘red¬ 
dish’ regions with enhanced [Sil]6717,6731/Ha line intensity 
ratio. Dark dust clouds are located to the west of the star 6r 
and to the east of the star Ir and around 5r, and also in the 
northern part of the region (these regions will be discussed 
in Section 4.3). 

The stratification of the emission is evident in the image 
showing the [Sll]6717,6731/Ha line intensity ratio (Fig. (91:), 
which is traditionally used as a shock indicator. At an aver¬ 
age electron density rie — 50—100 cm“^, as derived from the 
[Sii]6717/6731 flux ratio, /([Sii]6717,6731)//(Ha) > 0.4 is 
an evidence of a substantial contribution from the emissio n 
of gas behind the shock front (see, e.g., [Alien et al.ll2008l ). 
In Fig. [9]: regions with shock signatures are shown in red. 
(Note that the [Sll]6717,6731/Ha flux ratios may be under¬ 
estimated because the image of the region actually combines 
Ha and [Nil]6548,6584 lines (see Section 2), and this fur¬ 
ther supports our conclusion about the possible presence of 
shocks.) 

The most expressive shock manifestation is an X-shaped 
structure in the eastern part of the region between stars Ir 
and 5r (Fig. [9]:). However, high [Sll]/Ha intensity ratio in 


this area may also be due to low intensities of the lines re¬ 
sulting in large measurement errors. To separate the two fac¬ 
tors (shock excitation and observational errors), in Fig. [9]i 
we show a map of uncertainties of the derived [Sll]/Ha in¬ 
tensity ratio at each pixel. The comparison of Fig. [9]: and 
Fig. [fljl shows that the eastern part of the X-shaped struc¬ 
ture and the western part of the observed region are charac¬ 
terized by significant uncertainties of the [Sll]/Ha ratio and 
that the enhanced [Sil]/Ha ratio in these areas most likely 
arises due to very low intensity of the lines. Hence only the 
western part of the X-shaped structure can be considered as 
an indication of a shock front. 

The spectroscopic observations of the area confirm that 
gas emission in the vdB 130 region is actually a combi¬ 
nation of photoionized gas emission and (in some direc¬ 
tions) shocked gas emission. The [Nil]6548,6583/Ha line ra¬ 
tio along the slits (Fig. [TOl) exhibits high values in the same 
regions where the [Sli]/Ha line ratio is high. The average 
value of the [Nil]6548,6583/Ha line ratio is 0.40-0.45. The 
[Oiii]5007/H/3 ratio along the slits varies from 0.3 to 0.9. Lo¬ 
cal spikes of this ratio in Fig.[T0]may be associated with indi¬ 
vidual stars, e. g., 4r (slit Nl, -70...-150) or Hr and 12r (slit 
N4, 50-150) (see also Fig. [9^). The enhanced [Oiii]5007/H/3 
ratio in the region of the dust cloud to the east of the star Ir 
(slit N2, -70...70) may be due to large measurement errors 
of weak lines, as in the case of the [Sli]/Ha ratio. 

In Fig. im we relate our observed diagnostic dia¬ 
grams log ([OIIIJ/H/3) vs. log ([Nil]/Ha) and log ([OIIIj/H/3) 
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Figure 11. Diagnostic diagrams of the relations log ([OIIIJ/H/3) 
to log ([Nil]/Ha) (left) and to log ([SIIj/Ha) (right), observed 
with the spectrogra ms N1 and N3 (Fig . [OK)- Solid lines, accord¬ 
ing to the models of iKanffmann et all to03|l and iKewlev et alJ 
(|200lh and their oarametrization bv iKewlev et al.l (|2006l ). seoa- 
rate regions of pure photoionization excitation (below the black 
line), shock excitation (above both lines), and combined contri¬ 
bution of both mechanisms (between grey and black lines in the 
left panel). Typical measurement errors are shown in the bottom 
left corner of each panel. 


vs. log ([SII]/ H a) to theor etical photoioniza t ion m odels of 
iKewlev et al.l (|200lh and iKauffmann et al.l (|2003l ). Black 
lines mark the upper boundaries for possible line intensity 
ratios in the case of purely photoionization excitation mech¬ 
anism. A grey line in the left panel separates a domain with 
[Nil]/Ha and [Olll]/H/3 ratios typical of shocks (above the 
grey line) from the domain with the combined contribution 
of both mechanisms (between the black and grey lines). Ev¬ 
idently, in the studied region we observe mostly photoion¬ 
ization line excitation. However, a number of data points in 
the left panel fall inside the ‘composite’ domain. In Fig. [10] 
these points are marked by bigger circles on the profiles of 
the [Olll]/H/3 ratio. Comparing Fig.[T0](slits N1 and N3) and 
Fig. [9^, we see that domains with combined emission mech¬ 
anism coincide with the shell that is traced by the enhanced 
[Sll]/Ha ratio (Figs. [9l2,d), indicative of a possible contribu¬ 
tion from shocks. In particular, an analysis of spectroscopic 
observations confirms the combined emission mechanism in 
the western segment of the X-shaped structure, which is a 
part of the shell-like pattern visible in the [Sil]/Ha distribu¬ 
tion (see slit Nl, -80...-70 in Fig. (TOj). 

Thus, we conclude that an Hll region with possible 
shock signatures at its periphery is observed around a group 
of vdB 130 members. Unfortunately, we were able to perform 
a comprehensive analysis only for the Nl and N3 slits. The 
wavelength ranges for the other two slits include only [Oiil] 
and H/3 lines. 

We also used spectra Nl and N3 to estimate the gas 
metallicity. The [Oiil]4363 line that is needed for estimat¬ 
ing Te was too weak and we therefore could not apply the 
so-called direct method, which requires the knowledge of 
the electron temperature. An analysis of the gas metallicity 
based on the empirical method from IPilvugin &: MattssonI 
(|201lh calibrated by bright [Nil] and [Sli] lines yielded a 
close-to-solar oxygen abundance of 12 + log (0/H) = 8.6 zb 
0.15. 


4.3 Dust emission in the vdB 130 area 

Near-infrared images of the studied region in various Spitzer 
and Herschel bands are presented in Figs. [T]3,[7]3 andll2h- 
c. Infrared emission at 8 fim and 24 /im traces an unclosed 
shell with a diameter of 5 — 7 arcmin (about 3 pc) and a 
thickness from 0.5 to 1 arcmin (about 0.4 pc), which en¬ 
circles most of the stars of vdB 130. The shell is open in 
the direction of Cyg OBI association toward the east and 
partially overlaps with the molecular cloud toward the west 
(Fig. ( 7 ) 3 ). The brightest part of the shell is further referred 
to as Big Filament. Inside the shell there is an extended 
area of diffuse infrared emission with two blobs (hereinafter 
Blob E and Blob W), easily discernible on both IRAC and 
MIPS (24 jim) maps (Fig.[T2^,b,c). These blobs and the sur¬ 
rounding halo of bright infrared emission coincide in projec¬ 
tion with the head of the molecular cloud (Fig. [7 ]d) and are 
embedded into the region emitting at 160 /xm, shown with 
contours in Fig. 112b. Blobs W and E have sizes of ~ 1 arcmin 
(0.4 pc) and ~ 0.5 arcmin (0.2 pc), respectively. 

In Fig. 112b we show an expanded area of the 8 /xm 
emission map around the infrared blobs. Blob W (west¬ 
ern) is more extended and represents a nearly round patch 
of IR emission, with a sickle-like filament extending from 
its southern side. The fine structure with other filaments 
within the blob is visible at 8 /xm (as shown in more detail 
in Fig. [T2 ]d). The eastern blob (Blob E) is more compact 
and actually looks almost like a point source. These blobs 
and a sickle-like filament are also visible on a 24 /xm map 
(Fig. int)- 

Appearances of 8 and 24 /xm images are similar but not 
identical which is evident on the plot representing their ratio 
(Fig.[T2i). To produce this plot, we have convolved the 8 /xm 
map to the resolution of the 24 ^m map, u sing an appropri¬ 
ate kernel provided bv lAniano et 'al] (l201lli . While over most 
part of the studied region the ratio exceeds 1, there 

are few patches where it is smaller than 1, down to about 
0.2. In Fig.[T2]i these patches are encircled by white contours 
corresponding to F^jF^^ — 1. Smallest values of F^jF^^ are 
observed toward Blob E and Blob W, but there is another 
quite extended region where F^jF^^ is less than 1, located 
westward of Blob W, in the centre of a structure delineated 
by a curved part of Big Filament. 


A dissimilarity between 8 /xm and 24 /xm emission is 
a well-known featu r e of the ISM both i n our Galaxy 


(IWatson et ah I2OO8I: 

Deharvene et ahl I2OIOI) and in other 

galaxies ( 

jBendo et al. 

I2OO8I). In our Galaxy interiors of Hli 


regions tend to be brighter at 24 /xm while their borders are 
observed as ring-like structures on 8 /xm maps. It is tempt¬ 
ing to assume that dark spots in Figure 1121 1. where 8 /xm 
emission is suppressed relative to 24/xm emission, indicate 
the presence of Hli regions. Below we discuss whether or not 
other evidences may be presented in favour of this hypoth¬ 
esis. 

We start from the region encircled by Big Filament, 
which coincides with the region of ionized gas emission dis¬ 
cussed in the previous subsection. This region is relatively 
faint at 8 /xm and is filled with diffuse 24 /xm and Ha emis¬ 
sion (hereinafter we refer to it as Ha Blob as indicated in 
Fig. [12^). Big Filament itself is bright at 8/xm, less bright 
at 24 /xm, and dark in Ha emission. It is also clearly visible 
on 70 and 160 /xm images, obtained with Herschel^ and stays 
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Figure 12. (a) A map of the studied region at S/xm with overlaid contours of 500 /xm emission, (b) An expanded area of the S jiui map 
around the infrared blobs, (c) A map of the studied region at 24/xm with overlaid contours of 160 /xm emission, (d) A map of the ratio 
of fluxes at 8 p,m and 24 p,m. Dashed and solid lines indicate cuts for which flux ratio profiles are presented in Fig. 113k and Fig. 113b . 
respectively. 


discernible up to 500 jU m. The ring of infrared emission has 
been also mentioned in ISchneider et aP (|2007f ). 

The entire infrared emission pattern looks similar to 
Galactic Hli regions. To emphasize this similarity, in Fig¬ 
ure di we show how various infrared flux ratios change 
with the relative distance from the region centre along 
the cut shown in Figure fT^ with a dashed line. Also 
shown in Figure 113k are the same ratios in the Hil re¬ 
gion ROW 12 0, one of the best studied ‘infrared bubbles ’ 
in our Galaxy ([P eharven g et al.ll2009l : lAnderson et al.ir2Q10l : 
IPavlvuchenkov et al.l 1201^ 1. For all the three ratios trends 
are similar. The Fs/F 2 a ratio steadily grows from the region 
centre (marked by the smallest value of Fs/F 2 a) toward its 
border, while the T 24 /T 160 ratio decreases along the same 
cut down to 0.5 and then stays nearly constant. The ratio of 


Fs flux to T 160 flux exhibits the smallest variation staying 
approximately constant over the entire region. Both 8 /xm 
and 24 /xm maps have been convolved to a 160 /xm resolu¬ 
tion using kernels from lAniano et al.l (|201lh . 

Blob W is even more similar to RGW 120 in this respect 
as shown in Figure fT3b . with its interior being relatively 
brighter at 24 /xm than Ha Blob. Intensity ratios are shown 
for the cut indicated with a solid line in Fig. [T2jl. Blob E 
is not resolved well enough for a detailed comparison, but 
it is also faint at 8 /xm, bright at 24 /xm, and shows little 
variation in the Fq/Fiqq ratio. To ensure that this similarity 
is not bound exclusively to RGW 120, we have compared 
intensity ratios in our blobs with the cor responding ratios 
i n a f ew other infrared bubbles from the iGhurchwell et all 
(l200(t catalogue and found the same distribution on the 
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Figure 13. Profiles of intensity ratio for Ho; Blob (a) and Blob W 
(b). For comparison analogous profiles are also shown for a rep¬ 
resentative cut through the well-studied Hll region RCW 120. 

infrared fluxes. Thus, flux ratios in Spitzer and Herschel 
bands in all the three blobs shows spatial variations that 
are typical of the Hi I regions in the Galaxy. 

We have also performed an aperture photometry in 
Spitzer IRAC bands for several locations along Big Fila¬ 
ment, along the sickle-like feature in the southern part of 
Blob W, and also for Blob E and a nebulosity surrounding 
the star 4r (with the central source contribution subtracted). 
In each band surface brightness was measured for a circular 
region with a diameter of 5 arcsec. Background emission was 
estimated in extended regions with the lowest surface bright¬ 
ness within 100 arcsec from the studied objects. The regions 
listed above are indicated in Fig. [14^, while their infrared 
colours [3.6]-[4.5] and [4.5]-[5.8] are shown in Fig. [T^ . 

All the regions can be separated into three groups 
according to their relation to the domains described in 
I Ybarra et al.l (|2014l l. The first group comprises all the lo¬ 
cations in Big Filament along with the wl and w2 regions 
in the sickle-like feature. These regions have negative [3.6]- 
[4.5] colours an d [4.5]-[5.8] colours exceeding 3. According 
to the results of I Ybarra et al.l (|2014l T these colours are typ¬ 
ical of a photon dominated region (PDR) at a boundary of 


an Hil region with a moderate absorption (Ay ~ (4 — 5)"^) 
and log(G/nH) of the order of -1. Here G is the FU V ra¬ 
diation field in units of the Draine field ([Drain 3 I1978I T and 
nn is the gas density. Our estimate for the electron density 
of the order of 100 cm“^ implies G > 10. Note that in our 
picture all these regions are related to the boundary of the 
Hil region around vdB 130. 

The second group consists of Blob E and the w3^ w4, 
and w6 regions, that belong to Blob W. Their colours do not 
fall in to the domain considered in the work of I Ybarra et al.l 
(|2014I T with the [3.6]-[4.5] colour being positive and the 

[4.5] -[5.8 ] colour being sligh tly below 3. Extrapolating re¬ 
sults of I Ybarra et al.l (|2014l l , we may assume that absorp¬ 
tion in these directions is higher than 5"^ and log(G/nH) is 
lower than -2. We cannot expect that G is low in the region 
of bright near-IR emission, so lower log(G/nH) is proba¬ 
bly caused by higher density (consistent with higher absorp¬ 
tion). Note that extinction values inferred from I Ybarra et aU 
il2ni4h diagram agree with our estimates. 

Finally, the third group in our sample includes the w5 
region and the nebulosity around the star 4r. Location of 
the w5 region n early coincides with th e location of the star 
Ir. According to I Ybarra Ladal (|2009l L their colours ([3.6]- 

[4.5] ~ 0.6 and [4.5]-[5.8] ~ 2.2 — 2.4) can be explained by 
the emission of the low-temperature shocked H 2 gas. 

Thus, we may tentatively conclude that Ha Blob, 
Blob W, and Blob E and may be other regions with sup¬ 
pressed 8 /xm emission, also seen in Fig. [T2ji (like a region 
around the star 4r), are actually Hll regions surrounding 
young massive stars. In particular, the infrared emission 
properties of Ha Blob are entirely consistent with our con¬ 
clusion on the Hli region that surrounds the cluster vdB 130, 
based on the optical spectra. Blob W and Blob E can be ex¬ 
cited by the stars Ir (BIV) and 5r (B1V-2V), respectively 
(Fig. 12), while Ha Blob can be excited, among others, by 
stars Hr (BIV) and 12r, which are members of the cluster 
vdB 130. 

Being more extended and fainter in all infrared bands 
up to 160 /xm, Ha Blob is probably older than the other 
Hli regions in this area. We note that Big Filament which 
borders Ha Blob from the west, partially encircles a ring¬ 
like feature, seen on the map of the [Sll]/Ha ratio (Fig.[9l3), 
which is also visible as an absorption feature on the DSS 
image. The [Sil]/Ha ratio is somewhat enhanced along this 
ring-like feature relative to its interior, indicative of a shock 
(Fig. [3:). As we mentioned in the previous subsection, in the 
eastern part of the region we observe a clear signature of a 
shock. The sickle-like feature bordering Blob W from the 
west may have been formed by Ha Blob expansion, when 
the ionization front has started to engulf a parent clump for 
Blob W. 

The presence of Hli regions in the vicinity of vdB 130 
stresses the problem of its age overestimation. If Ha Blob 
is indeed excited by the cluster star(s), their interaction 
should have started quite recently, after the cluster itself had 
formed. It must be noted that the stars Ir and 5r, which we 
propose as exciting stars for Blob W and Blob E (as there 
are no other B stars in the blobs’ vicinity), may not be 
cluster members (see above). It is interesting that the three 
blobs are located at the tip of a long chain of millimetre 
sources that exten ds to the west far beyond Big Filament 
([Motte et al.ll2007l L This also hints that the CO cloud, also 
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Figure 14. (Left) Locations of the regions, for which the colour-colour diagram is presented, on the 8 p,m map. (Right) Colour-colour 
diagram for the regions shown on the left panel. The regions w3 and w4 and the regions E and w6 have been moved apart slightly to 
avoid overlapping. 


visible in the dust continuum, may be a remain of a bigger 
structure that has been destroyed by the expanding super¬ 
shell. It is also possible that the interaction between the 
pre-existing cloud and the expanding supershell has moved 
some molecular material in the vicinity of the cluster B- 
stars, initiating the formation of Ha Blob. 

We may ask then whether Blob W and Blob E forma¬ 
tion has been triggered by an expansion of Ha Blob. Their 
location hints at some causal link with the shell, but the 
available information is definitely not sufficient for any firm 
conclusion. 


5 CONCLUSIONS 

Detailed analysis of stellar populations and the interstellar 
medium in the vicinity of the cluster vdB 130 was performed, 
using optical observations taken with the 6-m telescope of 
SAO RAS and the 125-cm telescope of the Crimean Lab¬ 
oratory of SAI MSU as well as archival data from Spitzer^ 
Herschel Space Telescopes, 2MASS and UCAC4 projects. 

The following results are obtained. 

1. New revision of the stellar population of 
vdB 130 revealed 36 stars with proper motions dif¬ 
fering by less than 4masyr~^ from the centroid 
determined by iKharchenko et af] (|2013h ~ 

(—2.22,—4.39) masyr“^). The sample shows large scat¬ 
ter of colours (0.6 mag for faint stars) that can be 
considered as an evidence for substantial differential 
reddening across the cluster. The minimum colour 
excess E{B — V) ~ 0.79 zb 0.02 mag, and the appar¬ 
ent distance module {V — Mv)app ~ 15.0 zb 0.3 mag, 
whereas in NIR E{J — H) ^ 0.27 zb 0.02 mag and 
{K — MKs)app ~ 11.57 zb 0.3 mag. Since NIR data are 
corrupted by large and inhomogeneous extinction in much 
lesser extent then in optics, we consider cluster param¬ 
eters obtained from 2MASS photometric data as more 
reliable. Photometric data are consistent with young cluster 


age ranging from 5 to 10 Myrs and distance of about 
1.8 zb 0.3 kpc. 

2. The physical relation between vdB 130, which is lo¬ 

cated near the outer boundary of a supershell surround¬ 
ing the Cyg OBI association , and the cometary m olecular 
cloud with the size of 0.°4 ( Schneider et al is cor¬ 

roborated by (1) similar distances to the cluster and the 
association (1.5 - 1.8 kpc), which agree within the errors, 
and (2) the agreement of line-of-sight velocities of the cloud 
(Ulsr ~ — 1 -L 5kms“^), Cyg OBI stars, and hydrogen ion¬ 
ized by these stars (Ulsr ~ 3 — 9kms“^). The molecular 
cloud is located in the region influenced by the stellar wind 
and radiation of Cyg OBI and vdB 130 stars. 

3. Most vdB 130 stars are surrounded by an IR shell hav¬ 
ing a diameter of 3 pc, which is visible at all Spitzer bands 
(Figs.EJ) and[T2}. The inner part of the shell is partially filled 
with diffuse emission at 8 and 24 /xm, with 24 fim emission 
being brighter than 8 /xm emission. Such a relation between 
8 /xm and 24 /xm is typical for Galactic Hll regions, so we 
suggest that this shell marks the location of an ionization 
front in an Hil region, Ha Blob, which is also visible on our 
Ha and [Sll] maps and is probably produced by the hottest 
stars of vdB 130. At the eastern border of the Hll region 
encircled by the IR filament signatures of a weak shock are 
identified. 

4. There are at least two other Hi I regions in the area. 
Blob W and Blob E, which are more compact and mostly 
visible in infrared maps. Blob W and Blob E are excited 
by BIV and B1V-B2V stars, respectively. On the maps of 
the [Sll]/Ha ratio Blob W is delimited from the east by an 
arc-like feature which may be produced by the expansion of 
Ha Blob. 

5. Blob W and Blob E may represent the latest step in a 
star formation sequence in the area related to the pillar-like 
CO cloud, but in general the star formation history across 
the area seems to be quite complicated. 
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